clear
svmat e(b), names(col)
xpose, clear varname
rename v1 coef
gen i = _n
tempfile B
save `B', replace

clear
matrix V = vecdiag(e(V))
svmat V, names(col)
xpose, clear varname
rename v1 variance
gen i = _n
merge 1:1 i using `B', nogen

gen rel_time = .
forvalues i = 1/5 {
replace rel_time = -(6-`i') if i == `i'
}
forvalues i = 6/16 {
replace rel_time = `i'-6 if i == `i'
}

sort rel_time
gen ll = coef - invnormal(0.975)*sqrt(variance)
gen ul = coef + invnormal(0.975)*sqrt(variance)
keep rel_time coef ll ul
